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Abstract 

The cross section for the diffractive deep-inelastic scattering process ep — > eXp is mea- 
sured, with the leading final state proton detected in the HI Forward Proton Spectrometer. 
The data analysed cover the range xp < 0.1 in fractional proton longitudinal momen- 
tum loss, 0.08 < \t\ < 0.5 GeV~ 2 in squared four-momentum transfer at the proton 
vertex, 2 < Q 2 < 50 GeV 2 in photon virtuality and 0.004 < f3 = x/xp < 1, where 
x is the Bjorken scaling variable. For xp < 10~ 2 , the differential cross section has a de- 
pendence of approximately da/dt oc e 6 *, independently of xp, j3 and Q 2 within uncer- 
tainties. The cross section is also measured triple differentially in xp, (3 and Q 2 . The 
xp dependence is interpreted in terms of an effective pomeron trajectory with intercept 
ajp(0) = 1.114±0.018 (stat.)±0.012 (syst.) tomo (model) and a sub-leading exchange. 
The data are in good agreement with an HI measurement for which the event selection is 
based on a large gap in the rapidity distribution of the final state hadrons, after accounting 
for proton dissociation contributions in the latter. Within uncertainties, the dependence of 
the cross section on x and Q 2 can thus be factorised from the dependences on all studied 
variables which characterise the proton vertex, for both the pomeron and the sub-leading 
exchange. 
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1 Introduction 



Diffractive processes such as ep — > eXp have been studied extensively in deep-inelastic elec- 
tron^proton scattering (DIS) at the HERA collider [1-8], since understanding them in detail 
is fundamental to the development of quantum chromodynamics (QCD) at high parton densi- 
ties. The photon virtuality Q 2 supplies a hard scale for the application of perturbative QCD, so 
that diffractive DIS events can be viewed as processes in which the photon probes a net colour 
singlet combination of exchanged partons. A hard scattering QCD collinear factorisation the- 
orem [9] allows 'diffractive parton distribution functions' (DPDFs) to be defined, expressing 
proton parton probability distribution functions under the condition of a particular scattered 
proton four-momentum. The x and Q 2 dependences of diffractive DIS can thus be treated with 
a similar theoretical description to that applied to inclusive DIS, for example through the appli- 
cation of the DGLAP parton evolution equations [10]. 

Within Regge phenomenology, diffractive cross sections are described by the exchange of 
a leading pomeron (P) trajectory, as illustrated in figure 1. HI diffractive DIS data [3] have 
been interpreted in a combined framework, which applies the QCD factorisation theorem to 
the x and Q 2 dependences and uses a Regge inspired approach to express the dependence on 
the fraction xp of the incident proton longitudinal momentum carried by the colour singlet 
exchange. The data at low xp are well described in this framework and DPDFs and a pomeron 
trajectory intercept have been extracted. In order to describe the data at larger xjp, it is necessary 
to include a sub-leading exchange trajectory (M), with an intercept which is consistent [2] with 
the approximately degenerate trajectories associated with the p, to, a 2 and f 2 mesons. 



Figure 1: Schematic illustration of the diffractive DIS process ep — > eXp and the kinematic 
variables used for its description in a model in which the pomeron (IP) and a sub-leading (M) 
trajectory are exchanged. 

In many previous analyses, including [3], diffractive DIS events are selected on the basis 
of the presence of a large rapidity gap (LRG) between the leading proton and the remainder 

'For simplicity, the incident and scattered leptons are always referred to in the following as 'electrons', although 
the data studied here were obtained with both electron and positron beams. 
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X of the hadronic final state. A complementary way to study diffractive processes is by direct 
measurement of the outgoing proton using the HI Forward Proton Spectrometer (FPS) [7, 11] 
or its ZEUS counterpart [5]. Although the available statistics are smaller, the FPS method 
of studying diffraction has several advantages over the LRG method. In contrast to the LRG 
case, the squared four-momentum transfer at the proton vertex t can be reconstructed. The FPS 
method also selects events in which the proton scatters elastically, whereas the LRG method 
does not distinguish the elastic case from dissociation to excited systems Y with small masses 
My. The FPS also allows measurements up to higher values of xp than is possible with the 
LRG method, extending into regions where the sub-leading trajectory is the dominant exchange. 
Together, the FPS and LRG data thus provide a means of testing in detail the extent to which the 
variables xp, t and M Y associated with the proton vertex can be factorised from the variables 
(3 = x/xp and Q 2 describing the hard interaction. 

In this paper, a measurement of the cross section for the diffractive DIS process ep — > eXp 
using the FPS is reported. The t dependence is presented in the form of a differential cross sec- 
tion xp d 2 a/dt dxjp, from which the exponential slope of the t distribution is measured and its 
dependence on other variables studied. Diffractive reduced cross sections, a^ 4 \/3, Q 2 , xp, t) 
at \t\ = 0.25 GeV 2 , and a^ 3 \f3, Q 2 , xp) integrated over t, are also measured. These observ- 
ables are used to investigate the dependences on (3 and Q 2 , to extract the pomeron trajectory 
intercept from the xp dependence and to quantify the sub-leading exchange contribution. The 
data are also compared directly with the LRG measurement [3] in order to test the compatibility 
between the results obtained with the two measurement techniques and to quantify the proton 
dissociation contribution in the LRG data. 



2 Experimental Technique 

The data used in this analysis correspond to an integrated luminosity of 28.4 pb _1 and were 
collected with the HI detector in the years 1999 and 2000. In these years the HERA collider 
was operated at electron and proton beam energies of E e = 27.6 GeV and E p = 920 GeV, 
respectively, corresponding to an ep centre of mass energy of y/s = 319 GeV. 

2.1 HI detector 

A detailed description of the HI detector can be found elsewhere [12]. Here, the components 
most relevant for the present measurement are described briefly. 

Scattered electrons with polar angles 2 in the range 153° < 9' e < 177° are measured in 
a lead/ scintillating-fibre calorimeter, the SpaCal [13]. The energy resolution is a(E)/E 
7%/^ [GeV] © 1% and the energy scale uncertainty varies between 2.0% at a scattered elec- 
tron energy of E' e = 11 GeV and 0.5% at E' e = 27.6 GeV [14]. A Backward Drift Chamber 
(BDC) in front of the SpaCal is used to measure the electron polar angle with a precision of 

2 In the right-handed coordinate system used, the origin is at the nominal interaction point, with the +z axis and 
the polar angle 8 = in the direction of the outgoing proton beam (the 'forward' direction). The +x axis points 
towards the centre of HERA. Transverse momenta are measured with respect to the beam axis. 
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0.5 mrad and to suppress background where neutral particles fake the scattered electron signal. 
The SpaCal also has a hadronic section, with an energy scale known to a precision of 7%. 

The Central Tracking Detector (CTD), with a polar angle coverage of 20° < 9 < 160°, 
is used to reconstruct the interaction vertex and to measure the momentum of charged parti- 
cles from the curvature of their trajectories in the 1.15 T field provided by a superconducting 
solenoid. The finely segmented Liquid Argon (LAr) sampling calorimeter [15] surrounds the 
tracking system and covers the range in polar angle 4° < 9 < 154°. Its total depth varies with 
9 between 4.5 and 8 interaction lengths. The absolute hadronic energy scale is known with a 
precision of 4% for the measurements presented here. The hadronic final state is reconstructed 
using an energy flow algorithm which combines charged particles measured in the CTD with 
information from the SpaCal and LAr calorimeters [16]. 

The luminosity is determined with a precision of 1.5% by detecting photons from the Bethe- 
Heitler process ep — > epy in a crystal Cerenkov calorimeter, located at z = —103 m. 

The energy and scattering angle of the leading proton are obtained from track measurements 
in the FPS [7, 1 1]. Protons scattered through small angles are deflected by the proton beam-line 
magnets into a system of detectors placed within the proton beam pipe inside movable stations, 
known as Roman Pots. Each Roman Pot station contains four planes of five scintillating fibres, 
which together measure two orthogonal coordinates in the (x, y) plane. The stations used in 
this analysis approach the beam horizontally from outside the proton ring and are positioned at 
z = 64 m and z = 80 m. The detectors are sensitive to scattered protons which lose less than 
10% of their energy in the ep interaction and which are scattered through angles < 1 mrad. 

For each event, the leading proton energy and the proton scattering angles at the interaction 
point in the horizontal (x — z) and vertical (y — z) planes are obtained by applying transfer 
functions derived from the beam optics to the track parameters reconstructed in the FPS. The 
scattered proton energy is thus measured independently using the information in the horizontal 
and vertical planes. By comparison of these results, it is inferred that the energy resolution 
is around 6 GeV, independently of energy within the measured range, and that the absolute 
energy scale uncertainty is 0.5 GeV. The uncertainties in the reconstruction of the transverse 
momentum components p x and p y are quantified using a sample of elastic ep — > ep°p photo- 
production events with p° — > it + it~ decays. By comparing the FPS measurements with values 
reconstructed from the charged pions in the CTD, the resolution of the FPS is determined to be 
~ 40 MeV for p x and ~ 100 MeV for p y , dominated by the transverse momentum spread of the 
proton beam at the interaction point. The corresponding t-resolution varies over the measured 
range from 0.04 GeV 2 at \t\ = 0.08 GeV 2 to 0.08 GeV 2 at \t\ = 0.5 GeV 2 . The uncertain- 
ties in the transverse momentum measurements are 10 MeV for p x and 30 MeV for p y . The 
t-dependence measured in the FPS for the p° sample [17] is in good agreement with published 
HI data [18]. For a leading proton which passes through both FPS stations, the average overall 
track reconstruction efficiency is 20 ± 2%, corresponding to the product of the efficiencies in the 
individual scintillating fibre planes. The uncertainty on this efficiency is evaluated by varying 
the details of the reconstruction procedure, for example the number of fibres per plane which 
are required to register a track element. 
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2.2 Event selection and kinematic reconstruction 

The events used in this analysis are triggered on the basis of a coincidence between a track 
in the FPS, an electromagnetic cluster in the SpaCal calorimeter and a charged particle track 
providing an interaction vertex in the CTD. The trigger efficiency varies with the kinematic 
variables studied and is around 85% on average. 

Several selection criteria are applied to the data in order to suppress beam related back- 
grounds, background due to photoproduction processes and events in which the incoming elec- 
tron loses significant energy through QED radiation. The DIS selection criteria are summarised 
below. 

• The reconstructed z coordinate of the event vertex is required to lie within 35 cm 
(~ 3cr ) of the mean position. At least one track originating from the interaction vertex 
and reconstructed in the CTD is required to have a transverse momentum above 0.1 GeV. 

• The variables characterising the scattered electron, E' e and 9' e , are determined from the 
SpaCal cluster, linked to a reconstructed charged particle track in the BDC, and the inter- 
action vertex reconstructed in the CTD. The electron candidate is required to satisfy the 
criteria 155° < 9' e < 176.5° and E' e > 11 GeV. 

• The quantity E —p z , computed from the energies and longitudinal momenta of all recon- 
structed particles including the electron, is required to lie between 35 GeV and 70 GeV. 
Neglecting detector effects and QED radiation, this quantity is expected to be twice the 
electron beam energy for neutral current DIS events. 

The following requirements are applied to the leading proton measured in the FPS. 

• The measurement is restricted to the region where the FPS acceptance is high by requiring 
that the transverse momenta in the horizontal and vertical projections lie in the ranges 
—0.38 < p x < —0.24 GeV and \p y \ < 0.7 GeV, respectively, and that the fractional 
energy of the leading proton, E ' / E p , be greater than 0.9. 

• To suppress cases where a DIS event reconstructed in the central detector coincides with 
background in the FPS, for example due to an off-momentum beam proton (beam halo), 
the quantity E+p z , summed over all reconstructed particles including the leading proton, 
is required to be below 1900 GeV. Neglecting detector effects, this quantity is expected 
to be twice the proton beam energy for neutral current DIS events. 

The inclusive DIS kinematic variables, Q 2 , x and the inelasticity y, are reconstructed using 
the techniques introduced in [2]. In order to optimise the resolution throughout the measured 
y range, information is exploited from both the scattered electron and the hadronic final state 
according to 



y = vl + Vd(i - yd) ; Q 2 



tan 2 (^/2) 




(1) 
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Here, y e and y d denote the values of y obtained from the scattered electron only ('electron 
method') and from the angles of the electron and the hadronic final state ('double angle method'), 
respectively [19]. The analysis is restricted to the region 2 < Q 2 < 50 GeV 2 and 0.02 < y < 
0.6. 

With q, P and P' denoting the four-vectors of the exchanged virtual photon and the incom- 
ing and outgoing protons, respectively, further variables specific to diffractive DIS are defined 

as 

q-(P-P') Q 2 ... 

XP= q.p ' (3= 2q.(P-Pr () 

such that P can be interpreted as the fraction of the colourless exchange longitudinal momen- 
tum which is carried by the struck quark. Two different methods are used to reconstruct these 
variables. In the 'leading proton' method, xjp is reconstructed directly from the energy of the 
leading proton, such that 

X]P = l-E' p /E p ■ p= — . (3) 

X]p 

In the 'X-mass' method, the mass of the system X is first obtained from the hadrons recon- 
structed in the central detector using 

M 2 x = (E 2 -pl-pl-pl) had .y-, (4) 

Vh 

where the subscript 'had' represents a sum over all hadronic final state particles excluding the 
leading proton and yh is the value of y reconstructed using only the hadronic final state [20]. 
Including the factor yjyu leads to cancellations of many measurement inaccuracies. The diffrac- 
tive variables are then reconstructed using 

The results obtained with the leading proton and X-mass methods agree well in the low xjp 
range where both are applicable. The X-mass method is used for xjp < 0.006 and the leading 
proton method is used for xp > 0.006, the choice being made on the basis of which method 
provides the better resolution. 

The squared four-momentum transfer t — (P — P') 2 is reconstructed using the transverse 
momentum p t of the leading proton measured with the FPS and the best value of xjp as described 
above, such that 

/-* V l ■ / x2]pm P (to 

L (-min -. j 'ruin -. i \y) 

1 — Xjp 1 — X]p 

where |t m i n | is the minimum kinematically accessible value of \t\ and m p is the proton mass. In 
the analysis, the reconstructed \t\ is required to lie in the range 0.08 < \t\ < 0.5 GeV 2 . The 
final data sample contains about 3 300 events. 
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3 Monte Carlo simulation and corrections to the data 



Monte Carlo simulations are used to correct the data for the effects of detector acceptances and 
inefficiencies, migrations between measurement intervals due to finite resolutions and QED 
radiation. The reaction ep — > eXp is simulated using an implementation of the 'saturation' 
model [21] within the RAPGAP generator [22]. Following hadronisation using the Lund string 
model [23] as implemented in JETSET [24], the response of the HI detector is simulated in 
detail and the events are passed through the same analysis chain as is used for the data. Weights 
are applied to the generated events so that the important kinematic variable distributions are 
well described throughout the region of the measurement. 

The background from photoproduction processes, where the electron is scattered into the 
backward beampipe and a particle from the hadronic final state fakes the electron signature in 
the SpaCal, is estimated using the PHOJET Monte Carlo model [25]. This background is neg- 
ligible except at the highest y values and is 6% at most. The proton dissociation background, 
where the leading proton originates from the decay of a higher mass state, is estimated using an 
implementation in RAPGAP of the dissociation model originally developed for the DIFFVM 
Monte Carlo generator [26]. This background is negligible except at the highest x$> values, 
where it reaches 2.7%. Background also arises from random coincidences of DIS events caus- 
ing activity in the central detector with beam-halo protons giving a signal in the FPS. This 
contribution is estimated statistically by combining DIS events (without the requirement of a 
track in the FPS) with beam-halo protons from randomly triggered events. Subtractions of up 
to 7% are made as a function of the total reconstructed E + p z of the event. 

Cross sections are obtained at the Born level, using RAPGAP interfaced to the program 
HERACLES [27] to correct for QED radiative effects. The data are presented at fixed Q 2 , 
P, xjp and t values, with corrections applied for the influence of the finite bin sizes using a 
parameterisation of the '2006 DPDF Fit A to the HI LRG data [3] for the Q 2 , f3 and xjp 
dependences and the t dependences measured in this analysis at each xp value (see section 5.1). 



4 Systematic Uncertainties on the Measured Cross Sections 

Systematic uncertainties are considered from the following sources. 

• The uncertainties in the leading proton energy, its transverse momentum in the horizontal 
projection and that in the vertical projection are 0.5 GeV, 10 MeV and 30 MeV, respec- 
tively (see section 2.1). The corresponding average uncertainties on the and ov^ 
measurements are 5.7%, 6.0% and 3.3%. 

• The energy scale uncertainty of the SpaCal implies an error of between 0.5% and 2.0% 
(depending on the energy) on the E' e measurement, which leads to an average systematic 
error of 3.0% on the a® data points. Possible biases in the 9' e measurement at the level of 
±0.5 mrad lead to an average systematic error of 2.4%. 
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• The systematic uncertainties arising from the hadronic final state reconstruction are de- 
termined by varying the hadronic energy scales of the LAr and SpaCal calorimeters by 
4% and 7%, respectively, and the energy fraction carried by tracks by 3%. Each of these 
sources leads to an uncertainty in the a® measurements of typically 1.5%. 

• The model dependence of the acceptance and migration corrections is estimated by vary- 
ing the shapes of the distributions in the kinematic variables xp, (5 and t in the RAP- 
GAP simulation within the limits imposed by the present data. The xp distribution is 
reweighted by (l/xjp) ±ai , which leads to an average uncertainty of 2.0% in a® '. The 
P distribution is reweighted by f3 ±0A and (1 — (3) ±0A , leading to typical uncertainties of 
3.2%. Reweighting the t distribution by e ±t results in uncertainties of 2.5% on average. 

• The uncertainties related to the subtraction of backgrounds (see section 3) are at most 
2.7% for proton dissociation, 3.0% for photoproduction and 3.5% for the proton beam- 
halo contribution. 

• A 2.6% uncertainty is attributed to the trigger efficiencies (section 2.2), evaluated using 
independent triggers. 

• The uncertainty in the FPS track reconstruction efficiency results in an overall normali- 
sation uncertainty of 10% (see section 2.1). A further normalisation uncertainty of 1.5% 
arises from the luminosity measurement. 

• The extrapolation from the measured FPS range of 0.08 < \t\ < 0.5 GeV 2 to the region 
Kminl < \t\ < 1 GeV 2 covered by the LRG data [3] results in an additional systematic 
error of up to 5% for the data (see section 5.3). 



The systematic errors shown in the figures and tables are calculated as the quadratic sum of all 
contributions which vary from point to point, corresponding to average uncertainties of 12% for 
the (T, D( ^ data and 13% for a^ 3 \ The quoted errors do not include the overall normalisation 
uncertainty. 



5 Results and Discussion 
5.1 Cross section dependence on t 

The differential cross section d 2 a /dxpdt provides a measurement of the t dependence of 
diffractive DIS. This cross section is shown in figure 2a, multiplied by xp for convenience, 
for three values of t and six values of xp in the range xp < 0.1 and 0.08 < \t\ < 0.5 GeV 2 , 
integrated over 2 < Q 2 < 50 GeV 2 and 0.02 < y < 0.6. For each xp value, fits to the form 
xp d 2 a/dxpdt oc e Bt are shown in figure 2a. The extracted values of the slope parameter 
B are plotted as a function of xp in figure 2b and are listed in table 3. The HI results for B 
are consistent with ZEUS measurements [5], though the HI data are somewhat lower than the 
ZEUS data for x F < 0.02. 
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At low xp, the data are compatible with a constant slope parameter, B ~ 6 GeV 2 . In a 
Regge approach with a single linear exchanged trajectory, ap(t) = ap(0) + a'jpt, the slope 
parameter is expected to decrease logarithmically with increasing xp according to 

B = Bp — 2a'jp In xp , (7) 

an effect which is often referred to as 'shrinkage' of the diffractive peak. The degree of shrink- 
age depends on the slope of the pomeron trajectory, which is a'p ~ 0.25 GeV -2 for soft hadron- 
hadron scattering at high energies [28]. In contrast, vector meson measurements at HERA have 
resulted in smaller values of a'p, whether a hard scale is present [29, 30] or not [31]. Fits of 
the form of equation 7 are performed to the FPS data shown in figure 2b in the region where 
pomeron exchange is expected to dominate, namely to the three data points with 0.0009 < 
xp < 0.0094, for which the sub-leading exchange contribution is estimated to be at most 7% 
(see the fit results in section 5.2). A two parameter fit to the data in this range yields Bp = 
6.0±1.6 (stat.)+ 2 ;<S (syst.) GeV~ 2 and a'p = 0.02±0.14 (stat.)^ (syst.) GeV~ 2 . Extending 
the fit range to the interval 0.0009 < xp < 0.021, for which the contribution of the sub-leading 
exchange is at most 20% (section 5.2), results in Bp = 4.9 ± 1.2 (stat.)*^ (syst.) GeV -2 
and a'jp = 0.10 ± 0.10 (stat.)^ao? ( s y st GeV~ 2 . The data thus favour a small value of a', as 
expected in perturbative models of the pomeron [32]. However, the result a'p ~ 0.25 from soft 
interactions cannot be excluded. The results of these fits are summarised in table 1 . 



Range of Fit 


(GeV" 2 ) 


Bjp (GeV- 2 ) 


0.0009 < x F < 0.0094 
0.0009 < xjp < 0.021 


0.02 ± 0.0141°;^ 
0.10 ±0.010lg;$ 


6.0 ± 1.6+^ 
4.9 ± 1.2lJ;f 



Table 1 : The results of fits to the slope parameter data in two different ranges at low xjp in order 
to extract a'jp and Bp, together with their statistical (first error) and systematic (second error) 
uncertainties. 

A decrease of the slope B is observed towards the region of larger xjp > 0.03, where the 
contribution from the sub-leading exchange is expected to be significant (60% in the highest bin 
at xp = 0.076). This reduction of the slope parameter indicates that the size of the interaction 
region reduces as x F increases, reaching values of around 4 GeV -2 , characteristic of the spatial 
extent of the proton charge distribution. 

The t dependence of the cross section is also presented in figure 3 and table 4 in different 
regions of Q 2 and (3 for two xp intervals. No significant Q 2 or (3 dependence of the slope 
parameter B is observed for xp < 0.03. Within the uncertainties, the t dependence of the cross 
section in the pomeron dominated low x P region can therefore be factorised from the Q 2 and 
P dependences. Since there is also no strong evidence for any (3 or Q 2 dependence of B for 
xp > 0.03, the data are consistent with a similar factorisation for the sub-leading exchange 
contribution. 

5.2 Cross section dependence on x F and extraction of ckjp(O) 

The xp, (3 and Q 2 dependences of diffractive DIS are studied in terms of the diffractive reduced 
cross sections and Ur^. The former observable is related to the measured differential 
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cross section by [3] 



df3dQ 2 dx]pdt (3Q 4 



(l-Hy)-^ (4) (Ag 2 ,^t). (8) 



The reduced cross section is equal to the diffractive structure function F^ 4 \(3,Q 2 ,X]p,t) to 
good approximation in the relatively low y region covered by the current analysis, where the 
contribution from the longitudinal structure function F^ 4 \(3, Q 2 , x P , t) is small. Results for 

ov 3 ^ are obtained at a fixed value of \t\ = 0.25 GeV 2 , interpolating from the measured range 
0.08 < \t\ < 0.5 GeV 2 using the measured t dependence at each x P value (figure 2). Presenting 
the measurement at \t\ = 0.25 GeV 2 ensures that the systematic uncertainties associated with 
this interpolation are small. 

Figure 4 shows xpo®^ 1 for \t\ = 0.25 GeV 2 as a function of x P for different Q 2 and (5 
values (see also tables 5- 8). At medium and large /3 values, x P falls or is flat as a function 
of xp. Qualitatively this behaviour is consistent with a dominant pomeron contribution with an 
intercept a P (0) > 1 . However, x paf^ rises with x P at the highest x P for low f3 values, which 
can be explained by a contribution from a sub-leading exchange with an intercept a R (0) < 1. 



To describe the x P dependence quantitatively, a fit is performed to the structure function 
F^ 4 \ obtained by correcting for the small F®^ contribution using the results of the 
'2006 DPDF fit A' in [3]. A parameterisation of the form 

if (4) = fr(x P , t)Fjp(f3, Q 2 ) + n R ■ fjnixjp, t)F R (/3, Q 2 ) (9) 

is used. This parameterisation assumes a separate 'proton vertex' factorisation of the x P and t 
dependences from those on (3 and Q 2 for both the pomeron and a sub-leading exchange. The 
factors fjp and f R correspond to flux factors for the exchanges and are taken from the Regge- 
motivated functions, 

f F {x P ,t) = A P . 2Qf(t) _ 1 ; fm(x P ,t)=A R . , (10) 

Xjp Xjp 

assuming that the sub-leading exchange has a linear trajectory a R (t) = a R (0) + a' M t as for 
the pomeron. The values of A P and A p are chosen such that x P ■ J ™ m fjp,m,{x P , t) dt — 1 at 
xjp = 0.003, following the convention of [2]. The free parameters of the fit are the pomeron 
intercept «jp(0), normalisation coefficients F P (f3, Q 2 ) for the pomeron contribution at each of 
the nineteen (/3, Q 2 ) values considered, and a single parameter n p describing the normalisation 
of the sub-leading exchange contribution. 

A summary of the values assumed for the parameters which are fixed in the fits is given in 
table 2. The intercept «jr(0) of the sub-leading exchange is obtained from [2]. As in [2,3], 
the normalisation coefficients F M (f3, Q 2 ) for the sub-leading exchange in each (3 and Q 2 bin 
are taken from a parameterisation of the pion structure function [33]. The remaining fixed 
parameters describing the fluxes are taken from the present analysis. Averages of the two fits 
to the B{xp) data at low x P described in section 5.1 (table 1) are used to fix the pomeron 
parameters, B P = 5.5 GeV -2 and a'jp = 0.06 GeV" 2 . The behaviour of B(xjp) at large x P 
is sensitive to the parameters a' R and B R . Although the constraints are not strong, the data 
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are incompatible with the pair of values, a'^ = 0.9 GeV" 2 [34] and B M = 2.0 GeV" 2 [35], 
obtained from soft hadronic scattering data and applied previously in similar fits to F® data [2]. 
A good description of the slope parameter results over the full xjp range is obtained with B R = 
1.6 GeV~ 2 and a' R = 0.3 GeV -2 . A description based on this parameterisation is shown in 
figure 2b. 



Parameter 


Value 


a'jp 
Bp 
«k(0) 

a 'lR 

B]R 


0.06 GeV" 2 

5.5 ; 2 ;? GeV" 2 
0.50±0.10 

0.3 ±g;| GeV" 2 

1.6 GeV" 2 



Table 2: The values of the fixed parameters and their uncertainties, as used in the extraction of 
q;jp(0). Since a'jp and Bp are strongly anti-correlated when extracted from the data shown in 
figure 2b, they are varied simultaneously to obtain the errors on the fit results, as are a' M and 

B]R. 

The experimental systematic uncertainties on the free parameters are evaluated by repeating 
the fit after shifting the data points according to each individual uncertainty source described 
in section 4. A model dependence uncertainty is determined by varying the fixed parameters 
as described in table 2. The a'jp R and Bjpjp parameters are varied in the ranges given in the 
table, within which an acceptable description of the data is maintained, whilst requiring that 
a'jp and a' M lie between and the values describing soft hadronic scattering (0.25 GeV -2 and 
0.9 GeV~ 2 , respectively). The influence of neglecting the FjP^ contribution is also included in 
the model dependence uncertainty. 

As shown in figure 4, the fit provides a good description of the xjp dependence of the data 
(X 2 = 44 with statistical uncertainties for 51 degrees of freedom). Within uncertainties, the xjp 
dependence can therefore be factorised from the (3 and Q 2 dependences for each of the pomeron 
and the sub-leading contributions. 

The fit yields a pomeron intercept of 

a P (0) = 1.114 ±0.018 (stat.) ± 0.012 (syst.) (model) , 

the dominant uncertainty arising from the variations of a'jp and B F . This result for ajp(0) 
is compatible with that obtained from HI data measured using the LRG method [3] and with 
ZEUS measurements [5,6]. It is only slightly higher than the pomeron intercept describing soft 
hadronic scattering, ajp(0) ~ 1.08 [28]. However, if a'jp is set to the soft pomeron value of 
0.25 GeV -2 , ajp(0) increases to around 1.15. 

The result for the sub-leading exchange normalisation parameter is 

n R = [1.0 ± 0.2 (stat.) ± 0.1 (syst.) +J; 2 (model)] x 10~ 3 , 

the largest uncertainty arising from the variation of 0:^(0). The sub-leading exchange is im- 
portant at low (3 and high x F , contributing typically 60% of the cross section at the highest 
xtp = 0.08. 
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5.3 Cross section dependence on Q 2 and (3 

The reduced cross section a^ 3 \f3,Q 2 ,xjp) is defined as the integral of a^ 4 \/3, Q 2 , xjp, t) 
over the range |t m i n | < \t\ < 1 GeV 2 , which is the region covered by HI using the LRG 
method [3]. It is obtained here by extrapolating the FPS data from the measured range 0.08 < 
\t\ < 0.5 GeV 2 using the t dependence at each xjp value (section 5.1 and table 3). The extrapo- 
lation factor depends only weakly on xjp and is 1.7 on average, with an uncertainty of up to 5%. 
The measurement of xpo®^ is presented in figures 5-7 and tables 5- 8. The data are compared 
with predictions derived from the '2006 DPDF Fit A' to the LRG data presented in [3], with 
modifications as described in section 6.1. 

The Q 2 dependence of at fixed xjp and (3 (figure 5) is characterised by positive scaling 
violations (da?^/d\n Q 2 > 0) throughout the kinematic range, except possibly at the high- 
est (3 = 0.7. This observation is consistent with that from HI measurements using the LRG 
method [2, 3] and implies a large gluonic component to the DPDFs. As can be seen from the 
model comparison, the positive scaling violations may be attributed to the pomeron contribution 
even at the highest xjp values, where the sub-leading exchange is dominant. 

The dependence of oy 3 ^ on /3 is weak over most of the kinematic range (figure 6). Since the 
(3 dependence is determined in the quark-parton model by the diffractive quark densities, this 
implies that the quark densities do not decrease at the highest values of (3 studied. Indeed, cy 3 ^ 
clearly rises as (3 — > 1 at low Q 2 and xjp. Within the framework of DPDFs, this can be explained 
in terms of diffractive quark densities peaking at high fractional momenta at low Q 2 [2, 3]. The 
f3 dependence of diffractive DIS has also been interpreted in terms of the elastic scattering from 
the proton of colour dipoles produced by partonic fluctuations of the virtual photon [21,36,37]. 
In such models, the cross section at low and intermediate (3 values is dominated by qqg and 
qq fluctuations of transversely polarised photons, respectively. The rise of oy D( ' 3 ' ) as (3 — > 1 at 
low Q 2 has been interpreted in terms of qq fluctuations of longitudinally polarised photons [38], 
which are suppressed as Q 2 increases. 



6 Comparison with Other Measurements 
6.1 Comparison with HI large rapidity gap data. 

The FPS ar {3) data can be compared with HI measurements obtained using the LRG tech- 
nique [3], after taking into account the slightly different cross section definitions in the two 
cases. Firstly, the cross section ep — > eXY measured with the LRG data is defined to include 
proton dissociation to any system Y with a mass in the range M Y < 1.6 GeV, whereas Y is de- 
fined to be a proton in the cross section measured with the FPS. Secondly, if there are significant 
isospin-1 contributions to the sub-leading trajectory, charge-exchange reactions producing lead- 
ing neutrons are expected in the LRG measurement, which are not present in the proton-tagged 
FPS data. 

A point-by point comparison between the data obtained with the LRG and FPS meth- 
ods can be found in [3]. Here, the level of agreement is scrutinised in more detail in the range 
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xjp < 0.05, to which the LRG method is applicable. To make the comparison with a minimum 
of systematic uncertainty and to test for differences between the kinematic dependences of the 
two cross sections, the LRG measurement is repeated with an identical Q 2 , (3 and xp binning 
to that used for the FPS data. The ratio of the two measurements is formed for each (Q 2 , (3, 
xjp) point and the dependences of this ratio on each kinematic variable individually is studied 
by taking statistically weighted averages over the other two variables. Since the two data sets 
are statistically independent and the dominant sources of systematic error are very different, 
correlations between the uncertainties on the FPS and LRG data are neglected. 

The ratio of the LRG to the FPS cross section is plotted in figure 8 as a function of Q 2 , (3 
and xjp. The combined normalisation errors of 12.7% are not shown. Within the remaining 
uncertainties of typically 10% per data point, there is no significant dependence on (3, Q 2 or 
xjp. The ratio of overall normalisations, LRG / FPS, is 



the dominant uncertainties arising from the normalisations of the FPS and LRG data. This 
result is consistent with the prediction of 1.15^°;^ from the DIFFVM generator, where the total 
proton-elastic and proton dissociation cross sections are taken to be equal by default and their 
ratio is varied in the range 1:2 to 2:1 for the uncertainties [3,26]. 

Since the FPS measurement extends to larger x p values than the LRG measurement, the FPS 
data provide complementary constraints on the sub-leading exchange trajectory. The value of 
rip obtained in section 5.2 is compared with the similarly defined parameter obtained in [3] after 
dividing the latter by the factor 1 . 23 ± . 1 6 (equation 1 1 ) to account for the different M Y ranges 
of the two measurements. Since all other parameters describing the sub-leading trajectory are 
fixed to the same values in the two analyses, the ratio of rip results is then equivalent to the ratio 
of sub-leading exchange contributions in the two cross section measurements. The dominant 
model dependence uncertainties largely cancel when forming this ratio, which is 



where the first error is the combined statistical and experimental systematic uncertainty and the 
second is the residual model dependence uncertainty as defined in section 5.2. This result is 
consistent with unity, as expected for a dominantly isosinglet sub-leading trajectory {uo or /, 
rather than p or a exchanges). It is thus consistent with the conclusion from charge exchange 
cross section measurements obtained by tagging leading neutrons in DIS at x F = 0.1, which 
can be fully attributed to % exchange [7]. 

The predictions of the '2006 DPDF Fit A' to the HI LRG data [3] are compared with the FPS 
data in figures 5-7 after applying a factor of 1/1.39 (equation 12) to the sub-leading exchange 
contribution in the fit and an overall normalisation factor of 1/1.23 (equation 1 1) to account for 
the absence of the proton dissociation contribution in the FPS case. The FPS data are then well 
described in the region covered by the fit to the LRG data (Q 2 > 8.5 GeV 2 ). Extrapolating to 
lower Q 2 , the description remains reasonable. 

The good agreement, after accounting for proton dissociation, between the LRG and the 
FPS data confirms that the two measurement methods lead to compatible results, despite hav- 
ing very different systematics. The lack of any kinematic dependence of the ratio of the two 



a(M Y < 1.6 GeV) 
a(Y=p) 



1.23 ± 0.03 (stat.) ± 0.16 (syst.) , 



(ID 



o_n (LRG) 
on (FPS) 



1.39 ± 0.48 (exp.) ± 0.29 (model) , 



(12) 
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cross sections shows, within uncertainties, that proton dissociation with M Y < 1.6 GeV can be 
treated similarly to the elastic proton case. This supports the factorisation, for both the pomeron 
and the sub-leading exchange contributions, of processes occuring at the proton vertex from 
those describing the hard interaction, in terms of M Y as well as t (section 5.1) and xw (sec- 
tion 5.2). It also confirms that contributions from proton dissociation in the LRG measurement 
do not significantly alter the measured (3, Q 2 or x P dependences and hence cannot have a large 
influence on the diffractive gluon density or other information extracted from the LRG data. 



6.2 Comparison with ZEUS leading proton data. 

In figure 9 the FPS oy 3 ^ results are compared with those of the ZEUS collaboration, measured 
using their Leading Proton Spectrometer (LPS) [5] and also integrated over \t\ < 1 GeV 2 . 
The ZEUS data points are interpolated to the (3 and Q 2 values of this measurement using the 
dependences measured in [5]. There is very good agreement between the two data sets. The 
ratio of the ZEUS LPS to the HI FPS data averaged over the measured kinematic range is 
0.92 ± 0.04(stat.) ± 0.03(syst.) ± 0.15(norm.), which is consistent with unity taking into 
account the dominant normalisation uncertainties. Within the errors, there is no xp, (3 or Q 2 
dependence of the ratio. The sub-leading exchange contributes at a similar level at high x F and 
low (3 in both data sets. 



7 Summary 

A semi-inclusive cross section measurement is presented for the diffractive deep-inelastic scat- 
tering process ep — > eXp. The results are obtained using data taken with the HI detector at 
HERA, where the scattered proton carries at least 90% of the incoming proton momentum and 
is measured in the Forward Proton Spectrometer (FPS). The FPS data are in good agreement 
with those of the ZEUS collaboration obtained with their Leading Proton Spectrometer. 

The t-dependence is parameterised by an exponential function such that da/dt oc e Bt . The 
resulting values of the slope parameter B in the pomeron dominated range, xjp < 0.0094, are 
close to 6 GeV -2 and are independent of x$> in this range within errors, favouring an effective 
pomeron trajectory slope a'p which is close to zero. There is also no significant Q 2 or (3 depen- 
dence of B. The slope parameter decreases to around 4 GeV" 2 in the higher x F region, where 
an additional sub-leading exchange is found to contribute. 

The diffractive reduced cross section Ur^i^Q 2 ,xjp,t) is measured at \t\ = 0.25 GeV 2 . 
The xp dependence is described using a model which is motivated by Regge phenomenology, 
in which a leading pomeron and a sub-leading exchange contribute. The effective pomeron 
intercept describing the data is a P (<d) = 1.114 ± 0.018 (stat.) ± 0.012 (syst.) (model). 

The data are also analysed in terms of the diffractive reduced cross section a^ 3 \ obtained 
by integrating <Jr^ over the range |t min | < \t\ < 1 GeV 2 . At fixed x F , a relatively flat (3 
dependence is observed over most of the kinematic range. The data display scaling violations 
with positive da^ /<9 m Q 2 , except at the highest values of (3 ~ 0.7. 
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The FPS data are compared with the results of an HI measurement using events selected 
on the basis of a large rapidity gap (LRG) rather than a leading proton, which includes proton 
dissociation to states with masses M Y < 1.6 GeV. The ratio of the LRG to the FPS cross section 
is 1.23 ± 0.03 (stat.) ± 0.16 (syst.), independently of Q 2 , (3 and xjp within the uncertainties. 
Apart from this normalisation factor, the FPS and LRG measurements are in remarkably good 
agreement, despite having very different sources of systematic error. The magnitude of the 
sub-leading exchange component in the FPS data is compatible with that obtained from the 
LRG data, suggesting that charge exchange contributions in the latter are small. Within the 
present uncertainties, the HI diffractive DIS data are thus compatible with the factorisation of 
the variables x F , t and M Y associated with the proton vertex from the variables (3 and Q 2 , 
which describe the hard interaction, holding separately for the pomeron and for the sub-leading 
exchange trajectory. 

Acknowledgements 

We are grateful to the HERA machine group whose outstanding efforts have made this ex- 
periment possible. We thank the engineers and technicians for their work in constructing and 
maintaining the HI detector, our funding agencies for financial support, the DESY technical 
staff for continual assistance and the DESY directorate for support and for the hospitality which 
they extend to the non-DESY members of the collaboration. 

References 

[1] T. Ahmed et al. [HI Collaboration], Nucl. Phys. B 429 (1994) 477; 

T. Ahmed et al. [HI Collaboration], Phys. Lett. B 348 (1995) 681 [hep-ex/9503005]. 

[2] C. Adloff et al. [HI Collaboration], Z. Phys. C 76 (1997) 613 [hep-ex/9708016]. 

[3] HI Collaboration, "Measurement and QCD Analysis of the Diffractive Deep Inelastic 
Scattering Cross Section at HERA", DESY 06-049, submitted to Eur. Phys. J. C. 

[4] M. Derrick et al. [ZEUS Collaboration], Phys. Lett. B 315 (1993) 481; 

M. Derrick et al. [ZEUS Collaboration], Z. Phys. C 68 (1995) 569 [hep-ex/9505010]; 
J. Breitweg et al. [ZEUS Collaboration], Eur. Phys. J. C 6 (1999) 43 [hep-ex/9807010]; 
S. Chekanov etal. [ZEUS Collaboration], Eur. Phys. J. C 25 (2002) 169 [hep-ex/0203039]. 

[5] S. Chekanov et al. [ZEUS Collaboration], Eur. Phys. J. C 38 (2004) 43 [hep-ex/0408009]. 

[6] S. Chekanov etal. [ZEUS Collaboration], Nucl. Phys. B 713 (2005) 3 [hep-ex/0501060]. 

[7] C. Adloff et al. [HI Collaboration], Eur. Phys. J. C 6 (1999) 587 [hep-ex/981 1013]. 

[8] C. Adloff et al. [HI Collaboration], Nucl. Phys. B 619 (2001) 3 [hep-ex/0 106070]. 

[9] J. Collins, Phys. Rev. D 57 (1998) 3051 [Erratum-ibid. D 61 (2000) 019902] [hep- 
ph/9709499]. 

17 



[10] V. Gribov and L. Lipatov, Sov. J. Nucl. Phys. 15 (1972) 438 [Yad. Fiz. 15 (1972) 781]; 
V. Gribov and L. Lipatov, Sov. J. Nucl. Phys. 15 (1972) 675 [Yad. Fiz. 15 (1972) 1218]; 
Y. Dokshitzer, Sov. Phys. JETP 46 (1977) 641 [Zh. Eksp. Teor. Fiz. 73 (1977) 1216]; 
G. Altarelli and G. Parisi, Nucl. Phys. B 126 (1977) 298. 

[11] P. Van Esch et al, Nucl. Instrum. Meth. A 446 (2000) 409 [hep-ex/0001046]. 

[12] I. Abt et al. [HI Collaboration], Nucl. Instrum. Meth. A 386 (1997) 310; 
I. Abt et al. [HI Collaboration], Nucl. Instrum. Meth. A 386 (1997) 348. 

[13] R. Appuhn et al. [HI SPACAL Group], Nucl. Instrum. Meth. A 386 (1997) 397. 

[14] C. Adloff etal. [HI Collaboration], Eur. Phys. J. C 21 (2001) 33 [hep-ex/0012053]. 

[15] B. Andrieu et al. [HI Calorimeter Group], Nucl. Instrum. Meth. A 336 (1993) 499; 

B. Andrieu et al. [HI Calorimeter Group], Nucl. Instrum. Meth. A 350 (1994) 57. 

[16] C. Adloff et al. [HI Collaboration], Z. Phys. C 74 (1997) 221 [hep-ex/9702003]. 

[17] O. Karschnik, Ph.D. thesis, University of Hamburg, 2001, available from 
http://www-hl.desy.de/publications/theses-list.html. 

[18] S. Aid et al. [HI Collaboration], Nucl. Phys. B 463 (1996) 3 [hep-ex/9601004]. 

[19] S. Bentvelsen et al, Proceedings of the Workshop "Physics at HERA", eds. 
W. Buchmuller, G. Ingelman, DESY (1992), 23; 

C. Hoeger, ibid. 43. 

[20] A. Blondel and F. Jacquet, Proceedings of the Study of an ep Facility for Europe, ed. 
U. Amaldi, DESY 79/48 (1979) 391. 

[21] K. Golec-Biernat and M. Wusthoff, Phys. Rev. D 59 (1999) 014017 [hep-ph/9807513]; 
K. Golec-Biernat and M. Wusthoff, Phys. Rev. D 60 (1999) 1 14023 [hep-ph/9903358]. 

[22] H. Jung, Comput. Phys. Commun. 86 (1995) 147; 
(see also http://www.desy.de/~jung/rapgap.html). 

[23] B. Andersson, G. Gustafson, G. Ingelman and T. Sjostrand, Phys. Rept. 97 (1983) 31. 

[24] T. Sjostrand, Comput. Phys. Commun. 82 (1994) 74. 

[25] R. Engel and J. Ranft, Phys. Rev. D 54 (1996) 4244 [hep-ph/9509373]. 

[26] B. List and A. Mastroberardino, Proc. of the Workshop on Monte Carlo Generators for 
HERA Physics, eds. A. Doyle, G. Grindhammer, G. Ingelman, H. Jung, DESY-PROC- 
1999-02(1999) 396. 

[27] A. Kwiatkowski, H. Spiesberger and H. J. Mohring, Comput. Phys. Commun. 69 (1992) 
155. 



18 



[28] G. Jaroszkiewicz and P. Landshoff, Phys. Rev. D 10 (1974) 170; 

A. Donnachie and P. Landshoff, Phys. Lett. B 296 (1992) 227 [hep-ph/9209205]; 
J. Cudell, K. Kang and S. Kim, Phys. Lett. B 395 (1997) 311 [hep-ph/9601336]. 

[29] S. Chekanov et al. [ZEUS Collaboration], Eur. Phys. J. C 24 (2002) 345 [hep-ex/0201043]. 

[30] A. Aktas et al. [HI Collaboration], DESY 05-161, to appear in Eur. Phys. J. C [hep- 
ex/05 10016]. 

[31] J. Breitweg et al. [ZEUS Collaboration], Eur. Phys. J. C 14 (2000) 213 [hep-ex/9910038]. 

[32] J. Bartels and H. Kowalski, Eur. Phys. J. C 19 (2001) 693 [hep-ph/0010345]; 
I. Ivanov, N. Nikolaev and A. Savin, hep-ph/0501034. 

[33] J. Owens, Phys. Rev. D 30 (1984) 943. 

[34] W. Apel etal. [Serpukhov-CERN Collaboration], Nucl. Phys. B 154 (1979) 189 [Yad. Fiz. 
30(1979) 373]. 

[35] A. Kaidalov, Phys. Rept. 50 (1979) 157. 

[36] J. Bartels, J. Ellis, H. Kowalski and M. Wusthoff, Eur. Phys. J. C 7 (1999) 443 [hep- 
ph/9803497]. 

[37] J. Bartels, K. Golec-Biernat and H. Kowalski, Phys. Rev. D 66 (2002) 014001 [hep- 
ph/0203258]. 

[38] A. Hebecker and T. Teubner, Phys. Lett. B 498 (2001) 16 [hep-ph/00 10273]. 



19 



CM 
> 

<D 

o 

B 

S10 

•*-> 

~o 

X 

CM 

"D 

QL 



10 



: 0.0002 < x r < 0.002 


: 0.002 < x F < 0,006 


: 0.006 < x, < 0.014 






• H1 FPS 


• 

\ 




\ 

N 




— Fit~e Bt 




: 0.014 < x„ < 0.03 


'- 0.03 < x„ < 0.06 


0.06 <x,< 0.1 


• 

\ 


X 





0.2 0.4 0.2 0.4 0.2 0.4 

|t| (GeV 2 ) 



(a) 



^10 
0) 

O 

X 8 
00 



{ 
J- 



10 



-3 



• H1 FPS 

□ ZEUS LPS 
— Parameterisation 




10 



-2 



10 



MP 



(b) 



Figure 2: (a) The differential cross section xjp dV/dxjpdt measured in the kinematic range 
2 < Q 2 < 50 GeV 2 , 0.02 < y < 0.6 for different x%> intervals. The results of fits of the 
form xjpd 2 a/dxipdt oc e Bt are also shown, (b) The slope parameter B obtained from these fits, 
shown as a function of xp. The results obtained with the ZEUS LPS [5] and the parameterisa- 
tion of the HI data described in section 5.2 are also shown. The inner error bars represent the 
statistical errors and the outer error bars indicate the statistical and systematic errors added in 
quadrature. 
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Figure 3: (a,b) The differential cross section x$> dV/dxjpdt measured in different regions of 
(a) Q 2 and xp and (b) (3 andxp. The results of fits of the form xjpd 2 a/dxjpdt oc e Bt are shown. 
(c,d) The slope parameter B obtained from these fits, shown as a function of (c) Q 2 and (d) (3 
for two xp intervals. The inner error bars represent the statistical errors and the outer error bars 
indicate the statistical and systematic errors added in quadrature. 
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Figure 4: The diffractive reduced cross section x$>o r ' '((3, Q 2 , xjp, t), shown as a function 
of xjp for \t\ = 0.25 GeV 2 at different values of (3 and Q 2 . The inner error bars represent 
the statistical errors. The outer error bars indicate the statistical and systematic errors added 
in quadrature. An overall normalisation uncertainty of 10.1% is not shown. The solid curves 
represent the results of the phenomenological 'Regge' fit to the data, including both pomeron 
(IP) and sub-leading (FT) trajectory exchange, as described in section 5.2. The dashed curves 
represent the contribution from pomeron exchange alone according to the fit. 
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Figure 5: The diffractive reduced cross section xjp a r ([3,Q 2 , xjp) for \t\ < 1 GeV 2 , shown as 
a function of Q 2 for different values of xjp and (3. The inner error bars represent the statistical 
errors. The outer error bars indicate the statistical and systematic errors added in quadrature. An 
overall normalisation uncertainty of 10.1% is not shown. The solid curves represent the results 
of the 'HI 2006 DPDF Fit A' to LRG data [3], modified as described in section 6.1. The dashed 
curves represent the extrapolation of this prediction beyond the Q 2 range which is included in 
the fit. The dotted curves indicate the contribution of pomeron exchange alone in this model. 
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Figure 6: The diffractive reduced cross section xjp a r ([3,Q 2 , xjp) for \t\ < 1 GeV 2 , shown as 
a function of (3 for different values of xjp and Q 2 . The inner error bars represent the statistical 
errors. The outer error bars indicate the statistical and systematic errors added in quadrature. An 
overall normalisation uncertainty of 10.1% is not shown. The solid curves represent the results 
of the 'HI 2006 DPDF Fit A' to LRG data [3], modified as described in section 6.1. The dashed 
curves represent the extrapolation of this prediction beyond the Q 2 range which is included in 
the fit. The dotted curves indicate the contribution of pomeron exchange alone in this model. 
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Figure 7: The diffractive reduced cross section x jp ((3, Q 2 , xjp) for \t\ < 1 GeV 2 , shown as 
a function of xjp for different values of [3 and Q 2 . The inner error bars represent the statistical 
errors. The outer error bars indicate the statistical and systematic errors added in quadrature. An 
overall normalisation uncertainty of 10.1% is not shown. The solid curves represent the results 
of the 'HI 2006 DPDF Fit A' to LRG data [3], modified as described in section 6.1. The dashed 
curves represent the extrapolation of this prediction beyond the Q 2 range which is included in 
the fit. The dotted curves indicate the contribution of pomeron exchange alone in this model. 
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that for Y = p and \t\ < 1 GeV 2 , obtained from measurements using the LRG and FPS 
methods. The results are shown as a function of (a) Q 2 , (b) (3 and (c) xp, after averaging over 
the other variables. The lines represent the result of a fit to the data assuming no dependence on 
any of these variables. The inner error bars represent the statistical errors. The outer error bars 
indicate the statistical and systematic errors added in quadrature. Normalisation uncertainties 
of 12.7% are not shown. 
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Figure 9: The diffractive reduced cross section xpo r ((3, Q 2 , £p) for \t | < 1 GeV 2 , shown as 
a function of xp for different values of (3 and Q 2 . HI FPS data are compared with ZEUS LPS 
results [5]. The inner error bars represent the statistical errors. The outer error bars indicate 
the statistical and systematic errors added in quadrature. Normalisation uncertainties of around 
10% on each data set are not shown. 
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Table 3: The slope parameter B, extracted from fits to the data of the form da/dt oc e Bt in 
different regions of x P . The mean values of Q 2 , (3 and x P are also shown for each measurement. 
The first uncertainty given is statistical, the second systematic. 
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Table 4: The slope parameter B extracted from fits to the data of the form da/dt oc e Bt in 
different regions of xp, P and Q 2 . The mean values of these kinematic variables are also given 
for each measurement. The first uncertainty given is statistical, the second systematic. 
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Table 5: The diffractive reduced cross sections, x F a r ' measured at \t\ = 0.25 GeV 2 , and 
X]pa®^ integrated over |i m i n | < \t\ < 1 GeV 2 , measured at Q 2 = 2.7 GeV 2 and various (3 
and xjp values. The first uncertainty given is statistical, the second systematic. Normalisation 
uncertainties of 10.1% are not included. 
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0.0149 ± 0.0024+°-°^? 

— U.UU11 


5.3 


0.06 


0.0100 


0.0185 ± 0.00231°;°^ 


0.0138 ±0.0017l^| 


5.3 


0.06 


0.0220 


0.0240 ± 0.0031 + Hmq 


0.0173 ±0.0022+^028 

— U.UUzo 


5.3 


0.06 


0.0450 


0.0208 ± 0.0034+ l ;° l0 l2 ; 2 


0.0140 ± 0.00231^ 


5.3 


0.06 


0.0800 


0.0369 ± 0.0087 + Hnm 

— U.UUbU 


0.0246 ± 0.0057+°-°° 2 ? 

— U.UU41 


5.3 


0.15 


0.0011 


0.0260 ± 0.0029+^?! 

— U.UUlo 


0.0199 ±0.0022+Hm2 

— U.UUlO 


5.3 


0.15 


0.0040 


0.0202 ± 0.0021+Hmo 

— u.uuiz 


0.0153 ± 0.0016+Hmn 

— U.UU1U 


5.3 


0.15 


0.0100 


0.0222 ± 0.0030+^4 

— U.UU44 


0.0165 ± 0.0022+Hn^ 

— U.UUoo 


5.3 


0.15 


0.0220 


0.0243 ± 0.0034+°-°°^ 


0.0175 ±0.0024+°-°° 2 ° 

— u.uuzo 


5.3 


0.35 


0.0011 


0.0286 ± 0.0028^;^ 


0.0218 ±0.002ll°;°[! 2 ° 


5.3 


0.35 


0.0040 


0.0232 ± 0.0036 + S-^?I 


0.0177 ± 0.0027 + S$?i 

— U.UUlO 


5.3 


0.35 


0.0100 


0.0200 ± 0.0036l^^f 


0.0149 ±0.0027^;^ 


5.3 


0.70 


0.0011 


0.0460 ± 0.00371^^ 


0.0349 ± 0.00281^^? 


5.3 


0.70 


0.0040 


0.0419 ± 0.00561°;°^ 


0.0319 ± 0.0043+°;°°^ 



Table 6: The diffractive reduced cross sections, xjpa r measured at \t\ = 0.25 GeV 2 , and 
xjpa®^ integrated over |t min | < \t\ < 1 GeV 2 , measured at Q 2 = 5.3 GeV 2 and various (3 
and xjp values. The first uncertainty given is statistical, the second systematic. Normalisation 
uncertainties of 10.1% are not included. 
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Q 2 [GeV 2 ] 

L J 


P 




xjpa^ [GeV- 2 ] 

' L J 


D(3) 
Xjp<J r 


10.7 


0.02 


0.0100 


0.0325 ± 0.0054 + Hn64 


0.0243 ± 0.0040+^^ 


10.7 


0.02 


0.0220 


0.0345 ± 0.0040 + HnS 

U. UUJU 


0.0246 ± 0.00281°;^ 


10.7 


0.02 


0.0450 


0.0422 ± 0.00491^064 


0.0288 ± 0.0033lg;^ 


10.7 


0.02 


0.0800 


0.0502 ± 0.0075 + °-°n76 

— U.UU i D 


0.0330 ± 0.0050+° ^? 


10.7 


0.06 


0.0040 


0.0196 ± 0.00561°;°°^ 


0.0149 ± 0.0042+° °^ 

— U.UUU9 


10.7 


0.06 


0.0100 


0.0236 ± 0.0034+HnIs 


0.0177 ±0.0026 + Hnl 


10.7 


0.06 


0.0220 


0.0269 ± 0.0041+S-22?? 

— U.UU44 


0.0193 ±0.0029 + Hnl 

— U.UUoZ 


10.7 


0.06 


0.0450 


0.0329 ± 0.0054^;^! 


0.0224 ± 0.00371^;^ 


10.7 


0.06 


0.0800 


0.0278 ± 0.0087+°-°°^ 

— U.UUol 


0.0184 ± 0.0057 + ° ^4 

— U.UUo4 


10.7 


0.15 


0.0040 


0.0309 ± 0.0044+Hn9n 

U. UUj£U 


0.0232 ± 0.0033lg;^ 


10.7 


0.15 


0.0100 


0.0213 ± 0.0028 + Hnl4 

— U.UU44 


0.0160 ±0.0021 + HmI 

— U.UUo4 


10.7 


0.15 


0.0220 


0.0240 ± 0.0036 + Hnl 

— U.UUoo 


0.0173 ±0.0026 + Hn28 

— U.UUzo 


10.7 


0.15 


0.0450 


0.0254 ± 0.0047 + [! ^ 

— U.UUoO 


0.0174 ± 0.0032+^025 

— U.UUZO 


10.7 


0.35 


0.0011 


0.0382 ± 0.0055+^1 

— U.UUzl 


0.0292 ± 0.0042+°-°^ 

— U.UU1 1 


10.7 


0.35 


0.0040 


0.0292 ± 0.0051 + 2-{S?t 


0.0221 ± 0.0039 + ^^ 2 l 

— U.UU 14: 


10.7 


0.35 


0.0100 


0.0222 ± 0.00421^^3 


0.0166 ±0.0032 + S ^ 


10.7 


0.35 


0.0220 


0.0341 ± 0.0070+HnIq 


0.0246 ± 0.0051 + Hm6 

— U.UUOO 


10.7 


0.70 


0.0011 


0.0492 ± 0.00551°;°°^ 


0.0374 ±0.00411°;°^ 


10.7 


0.70 


0.0040 


0.0454 ± 0.0073lHoS 


0.0346 ± 0.0056lg;^ 


10.7 


0.70 


0.0100 


0.0339 ± 0.01021°;°°^ 


0.0254 ± 0.00771^053 



Table 7: The diffractive reduced cross sections, x F o r measured at \t\ = 0.25 GeV 2 , and 
xjpa®^ integrated over |t min | < \t\ < 1 GeV 2 , measured at Q 2 = 10.7 GeV 2 and various (3 
and xp values. The first uncertainty given is statistical, the second systematic. Normalisation 
uncertainties of 10.1% are not included. 
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Q 2 [GeV 2 ] 

L J 


P 




xjpa^ [GeV- 2 ] 

' L J 


D(3) 
Xjp<J r 


24.0 


0.02 


0.0220 


0.0425 ± 0.0088^;^ 


0.0306 ± 0.0063+^^ 


24.0 


0.02 


0.0450 


0.0497 ± 0.0090l°;°^ 2 


0.0341 ± 0.0062l°;™| 


24.0 


0.02 


0.0800 


0.0596 ± 0.0128+^086 


0.0394 ± 0.0085 + Ho58 


24.0 


0.06 


0.0100 


0.0264 ±0.0061+^ 


0.0200 ± 0.0046+°'°nl 

— U.UU4Z 


24.0 


0.06 


0.0220 


0.0386 ± 0.0070+°-°^ 

— U.UU45 


0.0276 ± 0.0050+° °^ 

— U.UUo4 


24.0 


0.06 


0.0450 


0.0334 ± 0.00751°,;°°$ 


0.0231 ± 0.0052+Hms 


24.0 


0.06 


0.0800 


0.0740 ± 0.0187l^;S 


0.0480 ±0.01241°;°°^ 


24.0 


0.15 


0.0040 


0.0252 ± 0.00681°,;°,°^ 


0.0190 ±0.00521^ 


24.0 


0.15 


0.0100 


0.0204 ± 0.0038 + °-°nIn 

— U.UU4U 


0.0152 ±0.0028 + Hn3n 

— U.UUoU 


24.0 


0.15 


0.0220 


0.0287 ±0.0051 + Hn!4 

— U.UU44 


0.0206 ± 0.0036 + Hn32 


24.0 


0.15 


0.0450 


0.0240 ± 0.0064+Hmq 


0.0164 ±0.0044 + °-°° 21 

— U.UUz i 


24.0 


0.15 


0.0800 


0.0414 ±0.0181 + SS 

— U.UU4o 


0.0274 ± 0.0120+^033 


24.0 


0.35 


0.0040 


0.0281 ± 0.0072+°-°°^ 


0.0214 ± 0.0055+2-^ 


24.0 


0.35 


0.0100 


0.0356 ± 0.0074+Hnfis 


0.0266 ± 0.0056 + Hntl 


24.0 


0.35 


0.0220 


0.0210 ± 0.0060+^^ 


0.0153 ±0.00431°,;°,°,^ 


24.0 


0.35 


0.0450 


0.0400 ± 0.0183lo;oo4g 


0.0274 ± 0.0126lo;oo34 


24.0 


0.70 


0.0011 


0.0535 ± 0.0125+S-^ 


0.0410 ± 0.0095 + Hn35 


24.0 


0.70 


0.0040 


0.0494 ±0.01061°;°,°^ 


0.0378 ±0.00811°;°°^ 


24.0 


0.70 


0.0100 


0.0238 ± 0.0071^Ho50 


0.0177 ±0.00541°,;°,°,^ 


24.0 


0.70 


0.0220 


0.0315 ± 0.0103^^053 


0.0231 ± 0.0074l°;™| 



Table 8: The diffractive reduced cross sections, x F o r measured at \t\ = 0.25 GeV 2 , and 
xjpa®^ integrated over |t min | < \t\ < 1 GeV 2 , measured at Q 2 = 24 GeV 2 and various (3 
and xjp values. The first uncertainty given is statistical, the second systematic. Normalisation 
uncertainties of 10.1% are not included. 
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